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Spontaneous entry of water molecules inside single-wall carbon nanotubes (SWCNT) 
has been confirmed by both simulations and experiments. Using molecular dynamics 
simulations, we have studied the thermodynamics of filling of a (6,6) carbon nanotube 
in a temperature range from 273 K to 353 K and with different strengths of the 
nanotube-water interaction. From explicit energy and entropy calculations using 
the two-phase thermodynamics (2PT) method, we have presented a thermodynamic 
understanding of the Filing behavior of a nanotube. We show that both the energy 
and the entropy of transfer decrease with increasing temperature. On the other hand, 
scaling down the attractive part of the carbon-oxygen interaction results in increased 
energy of transfer while the entropy of transfer increases slowly with decreasing the 
interaction strength. Our results indicate that both energy and entropy favor water 
entry in (6,6) SWCNTs. Our results are compared with those of several recent studies 
of water entry in carbon nanotubes. 
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I. INTRODUCTION 


Water has been a subject of interest for researchers for many decades. The ubiquitous 
presence of water makes it an interesting molecule to study in various physical conditions. 
Understanding the characteristics of water in conhned environment is important for design¬ 
ing various nanodevicesi*^ as well as for understanding various biological activities associated 
with different protein channels^>^. Both simulations^^- and experiments^^— have conhrmed 
spontaneous entry of water inside single-wall carbon nanotubes (SWCNTs). Inside the 
hydrophobic cavity of SWCNTs, water molecules exhibits many remarkable properties 
including highly ordered structuresi^il^, anisotropic rotational dynamics^ii^ and fast flow 
rates^. Recent studies have suggested potential applications of nanodevices based on carbon 
nanotubes for energy storage, desalination, voltage generation, flow sensing, fast flow devices 
etoi^^— . Inside narrow SWCNTs, water molecules are arranged in a solid-like structure to 
form single-hle chains of molecules with all dipoles pointing in the same direction along the 
nanotube axisl^ii^. This highly ordered one-dimensional structure of water molecules has 
been predicated to be thermodynamically stable^ for lengths up to millimeters at room 
temperature. Understanding the thermodynamics of water inside the hydrophobic cavity of 
a nanotube is interesting not only from the fundamental point of view but also to under¬ 
stand various biological activities where proteins undergo structural changes to control the 
presence of water inside various channels to perform specihc activities^i^^*^. 


In spite of many extensive studies^^^S.^— , the driving force of water entry into carbon 
nanotubes is not clear, mainly due to contradictory hndings from different studies. Using 
molecular dynamics (MD) simulations. Hummer et al^ hrst reported spontaneous entry of 
water molecules inside narrow SWCNTs and based on the energy distribution, proposed 
the driving force to be the gain in rotational entropy. Later Chandler et al— used a one¬ 
dimensional lattice-gas model to explore the thermodynamics of the hlling and emptying 
transitions of a hydrophobic tube. They explained the bimodal nature (empty and hlled) of 
the state of a nanotube immersed in water from an analogy with the liquid-vapor transition 
where it is energetically unfavorable to create a liquid-vapor interface but the entropy gain 
compensates for the energy increase. This simple model explains various observations from 
detailed atomistic MD simulations. However, due to its approximate nature, this model has 
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limited applicability. Later Vaitheeswaran et al.— computed energy and entropy of transfer 
using the temperature dependence of the occupancy probability for periodically replicated 
short segments of (6,6) SWCNT and concluded favorable entropy of transfer and unfavor¬ 
able energy of transfer for partial occupancy and the opposite behavior for fully occupied 
nanotubes. Recently, from energy and entropy calculations from atomistic MD simulation 
trajectories using the two-phase thermodynamics (2PT) method, we showed that water 
molecules inside a narrow SWCNT have higher rotational entropy as compared to the bulk 
and the gain in entropy is sufficient to compensate the increase of energy due to the loss of 
hydrogen bonding inside the cavity of the nanotube^^. Pascal et al— carried out a similar 
study for various nanotube diameters and found that the conhned water molecules have 
higher entropy up to a critical diameter of ~ 11 A and beyond this critical diameter, the 
energy of the conhned water molecules becomes favorable while the entropy decreases upon 
entry. Waghe et alM- computed the free energy of transfer from the occupancy probabil¬ 
ity for reduced carbon-oxygen interaction and extrapolated it to the actual carbon-oxygen 
interaction to argue that the entry of water molecules is driven by a favorable energy of 
transfer. Recently Garate et al— studied the free energy of transfer for various occupancies 
of SWCNT using thermodynamic integration and showed that for partial occupancy, the 
entry is driven by entropy, whereas the entry is energy driven for higher occupancies. 


It is important to understand the origin of such conhicting hndings and to develop a clear 
understanding of the thermodynamics of water entry inside the cavity of carbon nanotubes. 
In this report, we study the thermodynamics of single-hle water molecules conhned in a 
SWCNT at various temperatures and for diherent carbon-oxygen interaction parameters 
to develop a better understanding of the behavior. Based on these hndings, we try to 
understand the origin of the discrepancies among the results of earlier studies. 


The organization of the paper is as follows. In section HIl we give the simulation details 
and then present the results for the free-energy of transfer at diherent temperatures in 
section IIII Al In section IIIIBl we study the thermodynamics of conhned water for diherent 
carbon-oxygen interaction strengths. In section IIVI we compare our result with those of the 
other studies and explain the origin of the discrepancies among them. 
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II. METHOD OF SIMULATION 


To compute the free energy of filling, we have calculated the energy and entropy of water 
molecules in bulk and in conhnement inside a SWCNT. The differences between the energy 
and entropy of bulk and conhned water molecules give an estimate of the free energy of 
transfer. The energy of each water molecule can be directly computed from the simula¬ 
tion trajectories. To compute the entropy of water molecules, we have used the two-phase 
thermodynamics (2PT) method2^>^ based on the density of states which has been shown to 
give bulk water entropy values in excellent agreement with experiments^^. This method also 
has been applied to study the entropy of water and several organic liquids under various 
conditions^i^— . Details of the method to compute the entropy using this scheme have been 
presented in our previous study^. 

To study the thermodynamics of hlling at various temperatures, we have performed a 
series of MD simulations for a 54 A long (6,6) open-ended armchair SWCNT immersed in 
a bath of bulk water (~ 5000 water molecules) at temperatures ranging from 273 K to 353 
K. Interactions between various atom types were modeled using AMBER fflO force held—. 
Carbon atoms in the nanotube were modeled as Lennard-Jones (LJ) particles without any 
charge (ff03 atom type “CA”). We have used the TIP3P model — for water molecules and 
all OH bonds were constrained using the SHAKE algorithm. 


To systematically explore the effect of the carbon-water interaction on the thermodynam¬ 
ics of hlling, we performed the simulation with diherent carbon-oxygen interaction strengths. 
Following the procedure used by Waghe et al.—, we scaled the attractive part of the LJ in¬ 
teraction by a factor A such that: 


U(r, A) = 4e 




( 1 ) 


where e = eA^ and a' = a/A^/® are the modihed LJ parameters. We have studied the 
thermodynamics of hlling for four diherent values of A = 0.752, 0.785, 0.90,1.00. To make 
direct comparison with recent studies^!, we performed energy calculations for a 13.4 A long 
SWCNT solvated in bulk water (~ 3000 water molecules). Due to unfavorable interaction 
with the SWCNT for A < 1 values, water molecules do not stay long enough inside the 
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nanotube to produce sufficiently long trajectories to give converging velocity autocorrelation 
function data, making 2PT computations inaccurate. To overcome this problem, all the 
calculation for different A values were done for water molecules conhned inside a periodically 
replicated nanotube solvated in bulk water, which does not allow the water molecules to 
escape due to the presence of periodic images of the nanotube at both ends. In all cases, 
the nanotube was held hxed for the entire duration of simulation. 

Entropy and energy calculations were performed for occupancy = 5 which has been 
shown to be the minimum free-energy state for a 13.4 A long nanotube^I. To obtain the 
energy for different components of the system, we partition the potential energy into a sum 
over atoms^^. This is done by assigning equal portion of the energy to all participating 
atoms for vdW energy, bond energy, angle energy and dihedral energy. The per atom 
energy break up for the long-range Coulomb interaction part is accomplished by calculating 
the electrostatic potential at each charge site due to all other charges and multiplying the 
potential by the charge at that site. 

All the simulations were performed in the NPT ensemble at a constant pressure of 1 
atmosphere and temperatures ranging from 273 K to 353 K. Constant temperature and 
pressure were maintained using the Berendsen weak-coupling method as implemented in the 
PMEMD module of AMBER12^, with coupling constants of 2.0 ps and 2.0 ps, respectively, 
for temperature and pressure bath coupling. Long-range electrostatic interactions were 
computed using the particle-mesh Ewald summation scheme with a real-space cut-off of 
10.5 A. LJ interactions were computed using the same cut-off distance. The SWCNT was 
solvated with at least 15 A thick water shell in all 3 directions for all simulations except 
for the periodically replicated case in which there was no solvation shell in the 2 ; direction 
and both ends of the nanotube were connected across the boundary (see Figured]). 

After equilibration for 5 ns, coordinates and velocities of all water molecules were saved 
for 40 ps long trajectories with 4 fs saving frequency with integration time step of 1 fs. 
The density of states for entropy computation was calculated from the Fourier transform 
of the velocity autocorrelation function of water molecules. Translational and rotational 
components of the entropy were obtained by decomposing the total velocity into the center 
of mass velocity and the angular velocity. To get better statistical data, all calculations were 
performed for four independent trajectories and averaged values were obtained for different 
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physical quantities of interest. 


III. RESULTS 

Filling is governed by the difference between the chemical potentials of bulk and conhned 
water. The equilibrium between bulk water outside the nanotube and confined water inside 
can be altered by changing either of these chemical potentials. The chemical potential of 
bulk water can be changed by changing the temperature and the pressure or by adding salts. 
Changing the interaction strength between water and the carbon atoms that form the surface 
of the nanotube can also change the equilibrium between bulk and confined water. We have 
studied the effects of temperature and carbon-water interaction on the hlling behavior by 
computing the free-energy difference between confined and bulk water for different values of 
these quantities. 

A. Effects of Changing the Temperature 

We have computed the entropy and energy of conhned and bulk water molecules at differ¬ 
ent temperatures ranging from 273 K to 353 K. Figure [2] shows the variation of the transla¬ 
tional and the rotational entropy of transfer {AS = Sconf — Sbuik) at different temperatures. 
The translational entropy of conhned water is lower than that of bulk water molecules for 
all the temperatures considered here. In contrast, the rotational entropy of conhned water is 
always higher than that of bulk water over the temperature range considered here. However 
ASrot decreases with increasing temperature. This can be understood by considering hydro¬ 
gen bond breaking in bulk water at elevated temperatures. As the temperature increases, 
more hydrogen bonds break and bulk water molecules acquire relatively more rotational 
freedom and hence gain rotational entropy. On the other hand conhned water remains in 
the structurally same state with one hydrogen free and the other hydrogen bonded to one 
of the nearest water molecules in the chain. Hence, the diherence between the rotational 
entropies of conhned and bulk water molecules decreases with increasing temperature. 

The energy of transfer also exhibits a systematic variation with increasing temperature. 
Water molecules conhned inside the nanotube have lower energy and the energy of transfer 
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defined as AE = Econf — E^uik is a negative quantity. As the temperature is increased, the 
energy of transfer decreases (increases in magnitude), as shown in Figured This behavior 
can again be understood in terms of hydrogen bond breaking at elevated temperatures. 
In bulk water, hydrogen bonds can be broken more easily at higher temperatures due to 
increased kinetic energy. However water molecules inside the nanotube are less affected by 
thermal fluctuations due to strong conhnement and a smaller number of hydrogen bonds 
are broken as the temperature is increased. Upon increasing the temperature from 298 K 
to 333 K, the average number of hydrogen bonds per water molecule in the bulk changes 
from 3.55 to 2.68 while for the water molecules in conhnement, it changes from 1.80 to 1.56. 
Hence the change in energy of the conhned water molecules with increasing temperature is 
less as compared to that of the bulk water molecules. This causes the energy of transfer to 
decrease with increasing temperature. 

Combining the entropy of transfer with the energy of transfer, the free energy of transfer 
AA = AE — TAS is obtained. The temperature dependence of AA has been shown in 
Figure [3l This quantity is negative, consistent with the observation that the nanotube 
remains hlled during the simulation (see Table H]). It should be noted that these calculations 
were done for an open-ended nanotube immersed in a water bath which allows free exchange 
of conhned water molecules with those in the bulk. Hence the average occupancy of water 
molecules inside the nanotube changes slightly as the temperature is changed (see Table H]). 


B. Effects of Changing the Carbon-oxygen Interaction 

Filling of the carbon nanotube cavity by water molecules is extremely sensitive to the 
wall-water (carbon-oxygen) interaction strength. Hummer et ai— showed that reducing the 
depth of the carbon-oxygen LJ potential well by 0.05 kcal/mol results in the drying of 
the nanotube cavity for signihcant periods of time. We have studied the thermodynam¬ 
ics of Filing systematically as a function of the parameter A in Eq.([T]) which modihes the 
oxygen-carbon interaction strength. As shown in Figure 01 both translational and rotational 
entropies of conhned water increase with decreasing A, making the nanotube cavity entrop- 
ically more favorable. An increase in the entropy (both translational and rotational) with 
decreasing the strength of the attractive part of the interaction, as shown in Figure 0] and 
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Table m is physically reasonable. This observation is contrary to the results of Waghe et 
al.—, who found the entropy of transfer to be independent of A. 

The interaction energy of conhned water molecules for different A values has been shown 
in Table HIl As expected, scaling down the attractive part of the Lennard-Jones potential 
makes the nanotube cavity energetically less favorable and the energy of the water molecules 
inside the nanotube increases (become less negative) with decreasing A. The free energy of 
transfer shows the same trend, as shown in Figure [5l While the energy of transfer makes 
the cavity unfavorable for hlling, the entropy of transfer makes it favorable (see Figure H]). 
These two effects compete with each other to determine the free energy of transfer. The 
change in energy dominates over the change in entropy, and as shown in Table [ITl the free 
energy of transfer becomes more unfavorable with decreasing A. 


IV. DISCUSSION AND CONCLUSION 

When water molecules enter the nanotube, they lose hydrogen bonds, which increases 
their energy, whereas the vdW interaction with the carbon atoms in the nanotube decreases 
their energy. Reduced carbon-oxygen interaction makes the nanotube cavity unfavorable for 
water molecules and the nanotube remains empty for a considerable fraction of simulation 
time. In this case, a good estimate of the occupation probabilities, including the probability 
p(0) of having no water molecule inside the nanotube, can be made, which is not possible 
for the actual carbon-oxygen interaction. This allows the use of another method to compute 
the free energy of transfer for hlling. Logarithms of the ratio of the probabilities, p{N)/p{0), 
give estimates of the free energy of transfer AA^ for various occupancies. Waghe et al— 
computed the occupation probabilities p{N) from simulations of water inside a 13.4 A long 
nanotube of 8.1 A diameter with reduced carbon-oxygen interaction and extrapolated the 
results to the actual carbon-oxygen interaction. From these results, they concluded that it 
is energetically favorable, but entropically unfavorable for water molecules to be inside the 
nanotube, and hence the entry of water molecules inside a (6,6) SWCNT is energy driven. 
Recently Garate et al.—, using the thermodynamic integration method for a 31.9 A long 
(6,6) SWCNT, have shown that the entry of water is entropy driven for lower occupancies, 
but for the fully occupied state both energy and entropy favor conhnement. Favorable 


energy and entropy the fully occupied state of the nanotube is consistent with our result 
presented here, for both open ended (see Figures 2,3) and periodic nanotubes (see Figures 
4,5). However, these results are contrary to those of our previous study2^ where we found 
that the conhned water molecules have favorable entropy but unfavorable energy inside the 
nanotube. This difference lies in the method employed in computing the free energy of 
transfer. The method used by Waghe et ai— and Garate et alM‘ gives the total free energy 
of transfer of the whole system in going from zero occupancy to an occupancy state with N 
water molecules. In contrast the the free energy of transfer in our earlier work was computed 
as the difference in the energy and entropy of bulk and confined water molecules only and 
neglects the change of nanotube energy due to filling. 

For a 54 A long SWCNT at 300 K, the change in the nanotube energy due to Filing in 
the fully occupied state is ~ —2.04 kcal/mol per water molecule. The interaction energy per 
conhned water molecule inside a 54 A long nanotube at 298 K is —8.68 kcal/mol without 
including the change in the nanotube energy due to Filing. The interaction energy of water 
molecules “only” is higher than the bulk water molecules energy, —9.57 kcal/mol, at the 
same temperature. However, upon including the nanotube energy change due to Filing, the 
energy of a conhned water molecule becomes —10.72 kcal/mol indicating a favorable energy 
of transfer of —1.15 kcal/mol. 

This value qualitatively agrees with the values reported by Waghe et a7— (~-2.29 kcal/mol) 
and Garate et a/.— (~-1.21 kcal/mol) fTable HTB for open ended nanotube. The quantitative 
diherences might be due to diherent Filing states, diherent nanotube lengths or diherent 
water models used. The energy of conhned water molecules reported in this study includes 
the nanotube energy change per water molecule in all cases. The decomposition of the water 
molecules’ energy and nanotube energy due to Filing for all cases is given in the supporting 
information. We also observe that for partially occupied nanotubes, the energy of the con- 
Fned water molecules is higher than the energy in the bulk even after including the nanotube 
energy change. As the occupancy increases, the energy of the conhned water molecules be¬ 
comes more favorable, but the entropy of transfer becomes less favorable. This behavior is 
consistent with that observed by Garate et al. — who used the thermodynamics integra¬ 
tion method to compute the transfer free energy, energy and entropy. Results for diherent 
occupancies for a 13.4 A long periodic SWGNT are presented in the supporting information. 
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A comparison of the results for the energy, entropy and free energy of transfer obtained 
in different studies is shown in Table HU The results for AA obtained in different studies 
are similar, all indicating a favorable free energy of transfer, consistent with the observation 
that a nanotube immersed in a bath of water gets spontaneously hlled. The values of AE 
and AS obtained in this work are close to those reported by Garate et al. — . Waghe et al— 
also decomposed the free energy of transfer into entropic and enthalpic components. As¬ 
suming that the energy of transfer AEn and the entropy of transfer ASn in the expression 
A An = AEn —TASn are constant over the temperature range studied, a linear fit to AAn 
vs. T data was used to extract the values of AEn and ASn from the intercept and slope, 
respectively. The resulting values of AEn and ASn for reduced carbon-water interaction 
were extrapolated to A = 1 to obtain the values of these quantities for real carbon-water 
interaction. It was further assumed that the entropy of transfer is independent of carbon- 
water interaction and the energy of transfer was found to exhibit a linear dependence on the 
scaling factor A over the studied range. Based on these assumptions, Waghe et alM- argued 
that the confined water molecules have lower entropy as compared to bulk water molecules, 
in disagreement with the findings of this study and previous studies Results of our 

explicit calculations, presented in previous sections, show that the assumption of Waghe 
et al. about temperature independence of AE and AS in the temperature range studied 
does not hold true. This may be the reason for the discrepancy between the results and 
conclusions of Waghe et al. and those obtained in other studies^^i^i^. The reason for the 
differences between the results of Pascal et ali^ and those obtained in other studies remains 
to be explained.It should be noted, that the values of the energy and entropy of transfer are 
sensitive to the nanotube length and care should be taken to use same nanotube length while 
comparing the values presented here. For nanotubes of other diameters, both the energy and 
the entropy of transfer show systematic change as shown in the previous studies 

In summary, we have used fully atomistic MD simulations to study the effects of changing 
the temperature and the water-carbon interaction on the thermodynamics of the entry of 
water molecules in narrow carbon nanotubes. Both entropy and energy of transfer decrease 
with increasing temperature, keeping the free energy of transfer nearly constant. The en¬ 
tropy and the energy of water molecules inside the nanotube increase as we scale down the 
attractive part of the LJ interaction between the carbon atoms of the nanotube and the oxy- 
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gen atoms of water, with the energy increasing faster than the entropy. The change in the 
nanotube energy due to hlling is found to be important for understanding the differences be¬ 
tween the results of different studies of the driving forces for water entry into the nanotube. 
The results presented here will help in reconciling different views on the thermodynamics of 
water entry inside the hydrophobic channel of narrow carbon nanotubes. 
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FIG. 1. Simulated system. Left: 54 A long open ended single-wall carbon nanotube. Right: 
13.4 A long periodically replicated nanotube. Confined water molecules are shown in vdW repre¬ 
sentation. 



Temperature (K) 


FIG. 2. Variation of the translational, rotational and total entropies of transfer per water molecule 
(AS" = Sconf—Sbuik) with the temperature for a 54 A long open ended nanotube. Both translational 
and rotational entropies of transfer decrease with increasing temperature. 
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FIG. 3. Variation of the energy of transfer per water molecule {AE = E^onf — Ehuik) and the 
free energy of transfer per water molecule, AA, with the temperature for a 54 A long open ended 
nanotube. 



FIG. 4. Variation of the translational, rotational and total entropy of transfer per water molecule 
AS* as a function of A at T = 300 K for a 13.4 A long periodic nanotube with occupancy of 5 
water molecules. 
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FIG. 5. Variation of the energy and the free energy of transfer per water molecule, IS.E and 
A A = — FAS' respectively, as a function of A at T=300 K for a 13.4 A long periodic 

nanotube with occupancy of 5 water molecules. 


TABLES 


TABLE I. Energy, entropy and free energy of transfer per water molecule and average occupancy 
of the nanotube for different temperatures for a 54 A long open-ended nanotube. 


Temperature (K) AE (kcal/mol) TAS (kcal/mol) AA (kcal/mol.) Avg. Occupancy 


273.16 

-0.96 

± 

0.15 

0.73 

± 

0.11 

-1.69 

± 

0.26 

19.7 

285.15 

-1.01 

± 

0.15 

0.71 

± 

0.11 

-1.72 

± 

0.26 

19.4 

298.15 

-1.15 

± 

0.15 

0.58 

± 

0.14 

-1.73 

± 

0.29 

19.1 

313.15 

-1.28 

± 

0.15 

0.47 

± 

0.14 

-1.75 

± 

0.29 

18.9 

333.15 

-1.43 

± 

0.15 

0.35 

± 

0.14 

-1.77 

± 

0.29 

18.7 

353.15 

-1.48 

± 

0.16 

0.27 

± 

0.14 

-1.76 

± 

0.30 

18.3 
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TABLE II. Energy, entropy and free energy per confined water molecule inside a 13.4 A long 
periodic nanotube with occupancy of 5 water molecules for different values of the scaling factor A 
at T = 300 K. The values for bulk water are also shown for comparison. 


A 

E (kcal/mol) 

T Strans 

(kcal/mol) TSrot 

(k< 

;al/mol) E 

-TS 1 

(kc 

;al/mol) 

1.000 

-9.98 

± (0.16) 

3.88 

± (.20) 

1.69 

± 

(0.15) - 

-15.55 

± 

(0.51) 

0.900 

O 

O 

1 

± (0.16) 

3.93 

± (.20) 

1.74 

± 

(0.15) - 

-14.67 

± 

(0.51) 

0.785 

-8.42 

± (0.16) 

4.01 

± (.20) 

1.79 

± 

(0.15) - 

-14.22 

± 

(0.51) 

0.752 

1 

00 

± (0.16) 

4.05 

± (.20) 

1.83 

± 

(0.15) - 

-13.46 

± 

(0.51) 

BULK 

-9.45 

± (0.02) 

4.02 

± (.05) 

0.93 

± 

(0.01) - 

-14.40 

± 

(0.08) 


TABLE III. Comparison of the energy, entropy and free energy of transfer per water molecule for 
a (6,6) nanotube from various studies. The values are given for open ended nanotube only. 


Source AE (kcal/mol) TAS (kcal/mol) A A (kcal/mol) 


Pascal et al^ 

2.13 

3.06 

-0.93 

Garate et aL— 

-1.21 

0.43 

-1.63 

Waghe et alM- 

-2.29 

-0.76 

-1.53 

This work 

-1.15 

0.58 

-1.73 
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